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We study photon diffusion in a two-dimensional random packing of monodisperse disks as a 
simple model of granular media and wet foams. We assume that the intensity reflectance of disks 
is a constant r. We present an analytic expression for the transport-mean-free path I* in terms of 
the velocity of light in the disks and host medium, radius R and packing fraction of the disks, and 
the intensity reflectance. For the glass beads immersed in the air or water, we estimate transport- 
mean- free paths about half the experimental ones. For the air bubbles immersed in the water, l*/R 
is a linear function of 1/e, where e is the liquid volume fraction of the model wet foam. This throws 
new light on the empirical law of Vera et. al [Applied Optics 40, 4210 (2001)], and promotes more 
realistic models. 
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I. INTRODUCTION 

There is a good reason to study wave propagation 
in turbid or random media: Multiply scattered waves 
can probe temporal changes in physical systems p], [H, Q • 
Thus, light transport through fog 0], milky liquids, ne- 
matic liquid crysta ls H, g ranular media B,ui lp,[p , foams 
[Tol [Til [l2l [H, 113, Il5lll6j. and human tissue Il7| : propa- 
gation of elastic waves in the Earth's crust [iMUJl; acous- 
tic waves in the fluidized or sedimenting suspensions [2(| ; 
etc., have attracted much attention. 

In a turbid medium, light undergoes many scattering 
events before leaving the sample, and the transport of 
light energy is diffusive 0]. Therefore, the photon can 
be considered as a random walker. The transport-mean- 
free path I*, over which the photon direction becomes 
randomized, depends on the structural details of the 
opaque medium. Experimental techniques like diffuse- 
transmission spectroscopy (DTS) [2l| and diffusing-wave 
spectroscopy (DWS) [22| can be used to measure I*. In 
DTS, the average fraction T of incident light transmitted 
through a slab of thickness L is measured. The transport- 
mean-free path is then deduced from T oc I* / L. Utilizing 
the temporal intensity fluctuations in the speckle field of 
the multiply scattered light, DWS determines I* and the 
mean-squared displacement of the scattering sites due to 
time evolution, thermal motion, or flow. 

A plethora of light-scattering experiments show that 
light transport reaches its diffusive limit in granular me- 
dia [SSBH and foams [ISdEEBHEl, which 
means that photons perform a random walk. However, 
the mechanisms underlying this random walk are not elu- 
cidated. A wet foam is composed of spherical gas bubbles 
dispersed in liquid. A relatively dry foam consists of poly- 
hedral cells separated by thin liquid films. Three of them 
meet in the so-called Plateau borders which then define 
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tetrahedral vertices 23]. In their studies of foams with 
the liquid volume fraction e in the range 0.008 < e < 0.3, 
Vera, Saint- Jalmes, and Durian (llj observed the empir- 
ical law 

14 

J*«2iZ(— + 1.5), (1) 

£ 

where R is the average bubble radius. Recent studies 
of scattering from Plateau borders [ll|, H3 , vertices [25[ , 
and films [H H3, H, [H, H [M[ ; or transport effects such 
as total internal reflection of photons inside the Plateau 
borders [13, H3], have not yet clarified the empirical law 
of Vera et al.. For granular media, systematic measure- 
ments of the transport-mean- free path I* as a function 
of the refractive indices of grains and the host medium 
(air, water,...), grain size, and packing fraction, have not 
been performed. Menon et al. ((J determined I* ~ 15R 
for glass spheres of radius R = 47.5/im dispersed in air. 
For glass beads dispersed in water, Leutz et al. H found 
1* k UR ~ VoR for 80A«n < R sC 200/xm. Their samples 
had a packing fraction <j> k 0.64. Crassous H performed 
numerical simulations to find I* as a function of refrac- 
tive indices of the grain and host medium, but only for 
packing fraction (f> rj 0.64. 

It is instructive to consider simple or even toy models of 
granular media and wet foams, which allow an analytic 
access to the transport- mean- free path I*. Apparently, 
such models pave the way for deeper understanding of 
fascinating DWS experiments. In this paper, we consider 
two-dimensional packing of monosize disks. The disks 
are much larger than the wavelength of light, thus one 
can employ ray optics to follow a light beam or photon 
as it is reflected by the disks with a probability r called 
the intensity reflectance. We assume that the intensity 
reflectance is constant, and the velocity of light inside and 
outside the disks are c/rii n and c/n ou t, respectively. We 
show that the photon's random walk based on the above 
rules is a persistent random walk [26|, [H, [13] • Writing 
a master equation to describe the photon transport, we 
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find in Sec. If 11 Al the transport-mean-free path as 

r *R (^-i)(t^ + (¥) 2 ) (2) 

where and R denote the packing fraction and radius 
of disks, respectively. We further study our model by 
numerical simulation of the photon's random walk. We 
observe the overall agreement between our numerical and 
analytical estimates of the transport-mean-free path. 

For glass beads immersed in air or water, we find 
transport-mean-free paths about half the experimental 
ones |(| Q . For the air bubbles immersed in the water, 
we use Eq. © to derive I* as a function of the liquid 
volume fraction e = 1 — <j>. r m 0.20 is estimated as a 
weighted average of Fresnel's intensity reflectance. We 
find that in the range 0.08 < e < 0.15, our analytical re- 
sult agrees well with the relation I* w i?(0.11/e+2.37). In 
other words, we find that l*/R is a linear function of 1/e. 
Using the hybrid lattice gas model for two-dimensional 
foams and Fresnel's intensity reflectance, Sun and Hut- 
zler performed numerical simulation of photon transport 
and found I* w i?(0.26/e + 4.90) [H. Quite remark- 
ably, our analytic estimate of the transport-mean-frec 
path throws new light on the empirical law of Vera et 
al. and the numerical simulation of Sun et al.. 

Our article is organized as follows. In Section |TT] we 
introduce the two-dimensional packing of disks as a sim- 
ple model for a granular medium or a wet foam. Photon 
transport in a random packing of disks using constant 
intensity reflectance is discussed in Sec. IIIII Discussions, 
conclusions, and an outlook are presented in Sec. IIV1 

II. MODEL 

As a simple model for a two-dimensional disordered 
granular medium, wet foam, and bubbly liquid, we choose 
the random packing of circular disks. All non-overlapping 
disks have the same radius i?, and cover a fraction <j) 
of the plane. To address the photon transport in such 
medium, we have made the following assumptions: (i) 
Disks or grains, are much larger than the wavelength of 
light, thus one can employ ray optics to follow a light 
beam or photon as it is reflected by the disks with a 
probability r called the intensity reflectance, (ii) r is 
a constant, with no dependence on the incidence angle, 
(iii) Although disks of refractive index n;„ are immersed 
in a medium of refractive index n out , the incident and 
the transmitted rays have the same direction. In other 
words, we assume that the angle of refraction equals the 
angle of incidence, (iv) The velocity of light inside and 
outside the disks are c/ni n and c/n OMt , respectively. 

Our first assumption is inspired by the experiments 
[1, 0, B HI- Our second and third assumptions do not 
agree with Fresnel's formulas and SnelPs law, respec- 
tively. Consequently, our model does not consider to- 
tal internal reflection of rays. However, we deliberately 




c d 



FIG. 1: (a) Path of a photon moving in the host medium and 
hitting a disk with an incidence angle 7. (b) Path of a photon 
moving in a disk and hitting its surface with an incidence 
angle 7'. The step length inside the disk is 2Rcos"f' , where R 
is radius of the disk, (c) A photon impinging on a disk with 
an impact parameter s. Note s = J? sin 7. (d) The distance 
s' is related to the incidence angle 7' by s' = R sin 7'. 

adopt a step-by-step approach to photon transport in 
granular media, and will consider more realistic models 
later. 

As already mentioned, we model single photon paths 
in a packing of disks as a random walk with rules mo- 
tivated by ray optics, i.e., an incoming light beam is 
reflected from a disk surface with a probability r or it 
traverses the disk surface with a probability t = 1 — r. 
This naturally leads to a persistent random walk of the 
photons p6| . where the walker remembers its direction 
from the previous step [H, [34| . Persistent random walks 
are employed in biological problems (35| . turbulent dif- 
fusion |36[, polymers [37j |. Landauer diffusion coefficient 
for a one-dimensional solid [38[ , and in general transport 
mechanisms |39l. |40| | . More recent applications are re- 
viewed in [4l|. In the following section, we adopt the 
approach of |28l. [39| to study persistent random walk of 
the photons in a granular medium. 



III. PHOTON TRANSPORT IN A 
TWO-DIMENSIONAL PACKING OF DISKS 

A. Analytical treatment 

The photon random walk in a packing of disks con- 
sists of steps inside and outside the disks. We denote 
the average length of steps inside and outside the grains 
by Li„ and L out , respectively. We characterize each step 
by an angle relative to the £-axis. As Fig. [IJa) demon- 



3 



strates, on hitting a disk with an incidence angle 7, a 
photon moving in the host medium along the direction 
9 + tt + 2j will be either reflected to the direction 9, or en- 
ters the disk. The probability distribution of the random 
variable 7 (0 < 7 < it/2) is ^(7) = cos 7, see Appendix 
IA1 Similarly, a photon moving in a disk along the di- 
rection 9 + 7r + 27' and hitting its surface with an angle 
7', will be either reflected to the direction 9 or enters 
the host medium, see Fig. QJb). Quite remarkably, the 
probability distribution of the incidence angles 7 and 7' 
are the same, see Appendix [A] We can therefore con- 
clude that diffusion of the photons inside and outside the 
grains arc not inherently different. As will be shown in 
the following, we write a master equation to describe the 
photon diffusion inside (outside) the grains, utilizing the 
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step length L m (L out ) and velocity c/n m (c/n out ), and 
extract the diffusion constant Di n (D n ut)- According to 
the two-state model of Lennard- Jones [34( , the diffusion 
constant of photons in the granular medium is 



Djn fin^in f out D out? 



(3) 



where fi„ (font = 1 — fin) is the fraction of time that the 
photons spend inside (outside) the disks. 

We introduce the probability P n (x,y\9)dxdy that the 
photon after its nth step of length L along the direction 
9, arrives in the area dxdy at position x = (x,y). Then 
the following master equation expresses the evolution of 
P n (x,y\9): 



P n+ i(x,y\i 



P n (x - Lcos9,y - Lsm6\6 + tt + 2 7 ) F^dj + tP n (x - Lcos0, y - Lsin0|0). (4) 



r 



The first term on the right-hand side describes the re- 
flection of the photon with a probability r. The photon 
which has arrived at position (x — Lcos9,y — Lsui9) 
along the direction 9 + tt + 2j, changes its direction by 
an angle it + 27 according to the probability function 
-FX7) [42j. The second term describes the transmission 
with a probability t = 1 — r. The photon performs a bal- 
listic motion with step length L along direction 9 from 
position {x — L cos 9,y — L sin 9) to (x, y) . 

The diffusion constant follows from the evaluation of 
the second moment of P n {x, y\9) with respect to the spa- 
tial coordinates x and y. The probability distribution as 
an exact solution of the master equation Q is hard to 
obtain. However, there is a more direct method for the 
evaluation of the moments which employs the charac- 
teristic function P n (uj x ,ujy\m) associated with P n (x,y\9) 
S3: 



x kl y k2 P n (x,y\9)dxdyd9 



, z) k 1+ k 2 d kl+k2 Pn^\m = 0) 



,(5) 



j=0 



where k\ and k% are either zero or positive integers, 
P„(tf|m)=/ e lmB j f e iS x P n {x,y\9)dxdyd9, (6) 



and uj — (oj x ,uj v ). We are interested in the first and sec- 
ond moments of P n {x,y\9), thus we focus on the Taylor 
expansion 



P n (uj,a\m) « Q ,n(oc\m) + iu>LQ hn (a\m) 
^—Q2,n[a\m) + . . . , 



(7) 



where u> and a are the polar representation of the vector 
uj = (u> x ,ujy). From Eqs. ([5]) and it follows that 



{x)n 

(y)n 

(X 2 ) n 

(y 2 )n 



= ZQi,„(0|0), 

= LQi,„(||0), 

= £ 2 Q 2 ,„(0|0), 

= L 2 Q 2 ^\Q). 



(8) 



Fourier transforming Eq. we obtain 



P n+1 (u,a\m) = ike lka Jk{uL)V n {uj,a\k + m) 



k— — c 



x[r(-l) m+t F(2m + 2fc) + t], (9) 



where F(m) = \ fj 2 j2 e™ 7 F( 7 )d 7 , and 



J k (z) = 



e lzcos0 e' lke d9 



(10) 



is the fcth-order Bessel function. Since we are only inter- 
ested in the Taylor coefficients Qi >n (a\m) and Q2,n(o>\m), 
we insert Eq. ([7]) into Eq. ([9]). Using the Taylor expan- 
sion of the relevant Bessel functions Jk(z) (|fc| < 2) and 
Jfe(0) = (5 0j /c [iH and collecting all terms with the same 
power in uj, results in the following recursion relations for 
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the Q it n(a\m): 

Qo, n+ i{a\m) = [t + r(-l) m F(2m)]Q , n (a\m), 

Qi.n+Mm) = [t + r{-l) m F{2m)]Q 1 , n {a\m) 

+^Y~ [* + r{-l) m+1 F(2m + 2)] Q 0>n (a\m + 1) 

ia 

+— [* + 7-(-l) m_1 F(2m - 2)] Q ,n(a|m - 1), 
Q 2 ,n+i(a|m) = [t + r(-l) m F(2m)]Q 2 , n (a\m) 

+e- la [t + r(-l) m+1 F(2m + 2)]Qi,„(a|m + 1) 
+e M [i + r(-l) m - 1 F(2m - 2)]Q 1 ,„(a|m - 1) 

+\[t + r{-l) m F{2m)}Q , n {a\m) 

[t + r(-l) m+2 F(2m + 4)] Q^ n {a\m + 2) 

+— [* + r-(-l) m - 2 F(2m - 4)]Q 0>n (a|m - 2). 

(11) 



We solve this set of coupled linear difference equations 
using the method of the z-transform [34], Ej]. The 
z-transform Q(z) of a function Q n of a discrete vari- 
able n = 0,1,2,... is defined by Q(z) — Y^=oQnZ n ■ 
One then derives the z-transform of Q n +i simply as 
Q(z)/z — Q n —o/z. Note the similarities of this rule with 
the Laplace transform of the time derivative of a con- 
tinuous function [43| . The z-transform of equations JTT]) 
leads to a set of algebraic equations which immediately 
gives 



Qo(z,a\m) = 
Qi(z,a\m) = 



Qo,n=o(a\m) 



Q 2 {z,a\m) 



1 - z[t + r(-l) m F(2m)] ' 
Q\,n=o{ot\m) 



1 - z[t + r(-l) m F(2m)] 2(l - z[t + r(-l)" l F(2m)]) 

[t + r(~l) m + 1 F(2m + 2)] Q , n=0 (a\m + 1) e M [t + r(-l) m - 1 F(2m - 2)]Q Q , n=0 (q|m - 1) ' 
1 - z[t + r(-l) m + 1 F(2m + 2)] + 1 - z [t + r(-l)™- 1 F(2m - 2)] 

Qi,n=o{oi\m) zQ Qjn=0 {a\m) z 



l-z[t + r(-l)™F(2m)] 2 (l - z[< + r(-l)™F(2m)]) 2 1 - z[t + r(-l)"*F(2m)] 

[t + r(-l) m+1 F(2m + 2)] Qi(z, a|m + 1) + e M [t + r(-l) m-1 F(2Tn - 2)]Q x {z, a\m - 1) 

s 2M [t + r(-l) m F(2m - 4)] Q Q , n=0 (a\m - 2) 
4(1 - z[t + r(-l)™- 2 F(2m-4)]) 



e -2ia ^ + r (_i)m F ( 2m + 4)] Qo, n=0 {a\m + 2) 
4(1 - z[t + r(-l) m+2 F(2m + 4)]) 



(12) 



where F(m) = (1 — m 2 ) -1 cos(m7r/2), especially F(0) = 
1. The above expressions contain the sum of several 
terms whose inverse z-transform are readily accessible: 



n <-> 



na «-> 



1 - 


- z 




z 


(1 


-z) 2 ' 




1 


1 - 


- az ' 




2 


(1 


— az) 2 



(13) 



Here a is an arbitrary real number whose absolute mag- 
nitude is less than 1. 



For an arbitrary initial distribution Pq{x, y\9), the rel- 
evant function Qi, n (a\0) contains terms which are either 
constant or behave as a n with \a\ < 1. They are asso- 
ciated with the randomization of the initial distribution 
of the random walkers but are not essential for large n. 
According to Eq. 

{X) n = (y) n = (14) 

The behavior of the mean-square displacements is asso- 
ciated with Q2,m(a!|0), see Eq. ([HJ. We checked that for 
large n or in the long-time limit it is purely diffusive, i.e., 

(x 2 ) n = 2D x t, 

(y 2 ) n = 2D yT , (15) 
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where we introduced the time r = nL/v which passes 
when the random walker makes n steps at a speed v. 
We extract the diffusion constants from Q2n(0|0) and 
Ga,n(f|0): 



D x = D y = -lv{ 



2r 



1) 



(16) 



As already mentioned, we write the master equation 
(|4"|) to describe photon diffusion in the grains and in the 
host medium. Equation Q16p immediately leads to 



Di 



4 m n in { 2r 
4 nout 2r 



1), 



1). 



(17) 



The task is now expressing L in , L out , fi n and f out in 
terms of the model parameters R, <f>, r, and using the 
two-state model of Lennard- Jones to derive D m . First, 
we note that Lj„ =< 2i?cos7' >, where 7' is the inci- 
dence angle of photons moving in the disk, and here <> 
denotes averaging with respect to the probability distri- 
bution F'(-f') = COS7', see Appendix [Al and Fig. [ljb). 
Second, (f> = Li n /(Li n + L out ). Hence we find 



Lin 
Lout 



irR 

~2~' 
ttR 1 



(18) 



The evaluation of fc n is more exacting. Each step length 
Li n inside a disk takes a time r, n = Li n rii n /c. The prob- 
ability of m internal steps before leaving the disk is tr m . 
Hence the average time that a photon spends in the disk 
i s Em=o m ^ti' rm — Ti n r/t. The photon spends a time 
Tout = Loutnout/c before reaching a disk. Hence 



fir, 



n„r/t 



Tout +n n r/t nm -J— + n ou t - 

1-0 

TV,,,* Tlout j. 



Tout + T in r/t 



(19) 



Now we utilize Eq. ([3]) to derive the diffusion constant of 
photons in the granular medium as 



D, : 



ttRc (4? - !)( ih + i 1 ^-) 2 



(20) 



In two-dimensional space, the transport-mean-free path 
follows from 



(21) 



where c m is the transport velocity of light in the medium. 
To a first approximation 



(1 



(22) 



FIG. 2: Part of a packing of 10 4 disks, covering a fraction 
4> = 0.35 of the plane. 
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FIG. 3: The diffusion constant D m (in units of the disk radius 
R times the velocity of light c) as a function of the intensity 
reflectance r, for various packing fractions. Here fii n = 1.5 
and n ou t = 1.0 are assumed. Monte Carlo simulation results 
and D m (r) are denoted, respectively, by points and the line. 



Note that the velocity of light in the disks (the host 
medium) covering a fraction <f) (1 — 0) of the plane is 
cjfiin (c/ n ou t). From Eqs. (|20122p , we find the transport- 
mean- free path I* mentioned already in Eq. in the 
introduction. 



B. Numerical simulations 

We presented an analytic theory to calculate the diffu- 
sion constant of photons. Now we carefully examine this 
analytic result by performing numerical simulations. 

In order to generate random configurations of monodis- 
perse disks with a desired packing fraction, we compress a 
dilute system of rigid disks into a smaller space. Simula- 
tion methods based on a confining box generate a packing 
whose properties in the vicinity of walls differ from those 
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FIG. 4: The diffusion constant D m (in units of the disk radius 
R times the velocity of light c) as a function of the packing 
fraction <f>, for r = 0.1 and 0.4. Here m n = 1.5 and n out = 1.0 
are assumed. Monte Carlo simulation results and D m ((f>) are 
denoted, respectively, by points and the line. 



in the bulk. Hence we utilize the compaction method of 
Ref. 14511 which combines the contact dynamics algorithm 
[4r3.l47l| with the concept of the Andersen dynamics [4c| . 
This combined simulation method involves variable area 
of the simulation box with periodic boundary conditions 
in all directions. Due to the exclusion of side walls, the 
algorithm generates homogenous packings. 

We let photons perform a random walk in our packing 
of disks by applying the rules introduced in Sec. [TTJ For 
improving the speed of our ray tracing program, we adopt 
the cell index method commonly used in the molecular 
dynamics simulations (49j . The square simulation box is 
divided into a regular lattice of J x J cells. We maintain 
a list of disks in each of these cells. A photon moving 
in the cell j (1 ^ j ^ J 2 ) probably hits the disks in the 
cell j or its neighbor cells. Thus it is not necessary to 
check collision between the photon and all disks of the 
medium. 

Our computer program shrinks an initial dilute sample 
of 10 4 non-overlapping disks randomly distributed in a 
two-dimensional simulation box. In the course of shrink- 
ing the packing, the program saves snapshots of the grain 
positions if the packing fraction 4> & [0.15,0.25, ...,0.75], 
see Fig. [2] For each medium, the program takes 10 4 pho- 
tons at an initial position, and launches them in a direc- 
tion specified by angle 9q relative to the x-axis. Then it 
generates the trajectory of each photon following a stan- 
dard Monte Carlo procedure and evaluates the statistics 
of the photon cloud at times r G [7000, 7100, 9900] 
(in units of R/c). As detailed in Ref. [26|, we determine 
the diffusion constant D m from the temporal evolution 
of the average mean-square displacement of the photons: 
(x 2 + y 2 ) = 4D m T, 

For angles O € [20°, 40°, 320°], the simulation is re- 
peated for each intensity reflectance r € [0.1, 0.2, 0.9]. 
As a reasonable result, no dependence on the starting 
point and the starting direction is observed. In Fig. [3] 
we plot the diffusion constant D m (in units of the disk 
radius R times the velocity of light c) as a function of 
the intensity reflectance r, for the glass disks (rii n = 1.5) 



immersed in the air (n ou t = 1.0). For this medium, Fig. [4] 
shows the diffusion constant D m (in units of the disk ra- 
dius R times the velocity of light c) as a function of the 
packing fraction (f>, for the intensity reflectances r = 0.1 
and 0.4. The errorbars reflect the the standard deviations 
when we average over all diffusion constants for different 
starting positions and angles. We also performed simu- 
lations for the other examples (?ij n = 1.34, n ou t — 1.0), 
(n m = 1.0, nout = 1-34), (n in = 1.0, rw = 1.5), etc., 
but the results are not shown here. We observed the 
overall agreement between the numerical results and our 
theoretical value for D m . Quite remarkably, Eq. (j2"0|) in- 
volves no free parameters, but reasonably agrees with the 
numerical results in a wide range of 4>, r, m n , and n out - 



IV. DISCUSSION, CONCLUSIONS, AND 
OUTLOOK 



Diffusing-wave spectroscopy provides invaluable infor- 
mation about the static and dynamic properties ofgranu- 
lar media @, 0, H 1 and foams [U ED, El El Q M, El . 
The transport-mean-free path I* in terms of the micro- 
scopic structure, however, remains to be fully elucidated. 
In this paper, we consider a simple model for photon 
diffusion in a two-dimensional packing of disks. Our an- 
alytical result for I* provides new insights into the light 
transport, and promotes more realistic models. 

We have studied the photon's persistent random walk 
in a two-dimensional packing of monodisperse disks. We 
employed ray optics to follow a light beam or photon 
as it is reflected by the disks. We used a constant in- 
tensity reflectance r. Moreover we assumed that on hit- 
ting a disk, the incident and transmitted rays have the 
same direction. To achieve a better understanding of 
photon diffusion in granular medium and wet foams, we 
are extending our studies by considering Fresnel's formu- 
lae for the intensity reflectance, Snell's law of the refrac- 
tion, and the three-dimensional packing of polydisperse 
spheres. We are also improving our estimate of I* by tak- 
ing into account the distribution of photons' step length, 
and the transport velocity of photons [5(j. We discuss 
these points in the following. 

Many two-dimensional systems offer a rich and un- 
expected behavior. For an experimental observation of 
photon diffusion in a two-dimensional packing, a set of 
parallel fibres can be used. Photons injected in a plane 
perpendicular to the axis of fibres perform a planar dif- 
fusion. One can also drill parallel cylinders in a host 
medium, and fill all the cylinders with a liquid: I* de- 
pendence on the refractive index m n can be measured. 
Note that a ray maintains its polarization state on hit- 
ting a disk, and Fresnel's intensity reflectance depends on 
the polarization state: Quite interesting, the transport- 
mean-free paths for the transverse electric and transverse 
magnetic polarizations, are different. Inspired by the rich 
optics of a two-dimensional packing and following a step- 
by-step approach to a real system, our attention is natu- 
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FIG. 5: The probability distribution G(L ou t) as a function of 
Lout/R for various packing fraction <j>. Inset: The same plot 
in the logarithmic scale. Note that G(L ut) decays exponen- 
tially. 



rally directed to estimate I* of a i/iree-dimensional gran- 
ular medium or wet foam. We speculate that one can 
multiply Eq. ^ with a factor of about 1 to estimate the 
transport-mean-free path for a three-dimensional pack- 
ing. Apparently, only a detailed study will approve or 
disapprove our speculation which is based on the follow- 
ing observation: To understand the role of liquid films 
for light transport in dry foams, we studied the two- and 
three-dimensional Voronoi foams [13, Hl| ■ The interest- 
ing result is that the transport-mean-free path for these 
dry foams are determined by the same factor (1 — r)/r 
for a constant intensity reflectance r in spite of the differ- 
ence in dimension: Voronoi( r ) ~ l-10i?(l — r)/r and 
^3D Voronoi ( r ) ~ l-26i?(l — r)/r, where R denotes the av- 
erage cell radius [Ell ]. 

Our "mean- field" theory presented in Sec. IIIIA1 relics 
on the average step lengths Li n and L out . We note that 
Li n = 2i?cos7' and the probability distribution of the 
incidence angle 7' is F'(j') = cos 7', see Appendix [Al and 
Fig- Gib). It follows that the probability distribution of 
Lin is G(L m ) = L in /(2Ry/4R* - L? n ), L m = wR/2, and 
L?„ = / Lf n G(Un) dL ln = 8R 2 /3. Fig. [5] delineates the 
probability distribution G{L out ) as a function of L out /R 
for various packing fraction (f>. After reaching its pro- 
nounced maximum, the distribution G{L out ) decays ex- 
ponentially. The pioneering work of Heiderich et al. [52| 
suggests that the broadness of G(Li n ) and G(L OMt ) af- 
fects the value of I* by a factor about 1. 

In two-dimensional space, I* and D m are related as 
D m = l*c m /2, where c m is the transport velocity of light. 
In a medium composed of dielectric spheres comparable 
to the light wavelength A (Mie scatterers), the transport 
velocity differs by an order of magnitude from the phase 
velocity [53[. When spheres are much larger than the 
light wavelength, the difference between the two veloc- 



ities becomes unimportant. However, our first approxi- 
mation to c m [Eq. (|22[) ] can be improved by considering 
the real (rather than infinite) size parameter of disks, and 
the pair-correlation function of disks. 

We assumed that on hitting a disk, the incident and 
transmitted rays have the same direction. Note that in 
the case n out > rii n , a photon moving in the host medium 
and hitting the disks, experiences an average scattering 

angle J^ 2 {^ ~ %l)F{l)dl = 2 (in radians) due to the 
reflections. Taking into account inequality of the inci- 
dence and refraction angles, and the total internal re- 
flection of rays with incidence angle greater than "f c = 
arcsin(rij n /n ou t), the average scattering angle due to the 
transmissions is J Q 7c [arcsin(n out /m n sin 7) — 7^(7)^7 — 
0.132 (in radians), where n out = 1.34 and rij n = 1 are 
assumed. Reflections are more efficient than transmis- 
sions in randomizing the direction of photons. Thus the 
assumption that the incident and transmitted rays have 
the same direction, leads to a plausible estimate of I*. 

Fresnel's intensity reflectance depends on the incidence 
angle and the polarization state of the light. However an 
average (constant) intensity reflectance r may describe 
the photon diffusion. To estimate r, we first consider 
the case n ou t > rii n and define the critical angle 7 C = 
arcsin(nj n /n ou t). The probability distribution of the in- 
cidence angles 7 and 7' are F(j) = cos 7 and F'(j') — 

cos 7'. We define r out ^ in = J*' 2 r out ^ in (^)F(^)d^ . 
Here out — > in indicates a ray incident from the host 
medium onto the disk. The incident electric field can 
be decomposed into a component parallel (p) to the 
plane of incidence, and a component perpendicular (s) 
to this plane. To address the photon diffusion in a three- 
dimensional random packing, as reflectance r out ^i n (j) 
we have taken an average over the p and s polarizations: 
Tout-nnil) = 0-5[ r out->in (l, p) + r ut^in{l 1 s)] . Note that 
for 7 > 7c, the Fresnel's intensity reflectances are 1. 

Similarly, we define r m ^ out = J^ /2 0.5[r ln ^ out (j' ,p) + 
rin^outh' , s)]F' ('j')d-y' , and then use r = 0.5(r out ^ in + 
rin^out) in Eq. © to estimate I*. 

In the case n ou t < rii n , the appropriate critical angle 
is 7 C = arcsin(n ou t/rii rl ). Here < 7 < n/2 and 
< 7' < 7 C . The fact that 7' < 7 C ensures that a 
photon moving in the disk is able to enter the host 
medium, since the Fresnel's intensity reflectance is 1 for 
7' > 7c- An extension of Appendix [A] and numerical 
simulations indicate that the probability distribution 
of the incidence angles 7 and 7' are F(j) = cos 7 and 
F'(y) = cos 7'/ sin 7c, respectively. Wc define 

Jq /2 0.5[r out ^ in (j,p) + r out ^ in (j,s)}F(j)d'f and 
rin^out = Jq C 0.5[r ilwotlt (7',p)+ri„^ otlt (7',s)]F'(7')d7' 
to estimate r = 0.5(r otlt ^ n + ri n ^ out ). In the limit 
riin = n out we find r = 0. In this limit the photon 
transport is ballistic and Eq. ^ correctly predicts 

r->oo. 

For the glass disks (rii n = 1.5) immersed in the air 
(n-out = 1-0), we estimate r ~ 0.12. Thus for a packing 
fraction <fr = 0.64, we find I* w 7R. For the glass disks 
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Analytical Estimate 
0.11/E+2.37 
0.26/E+4.9 
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FIG. 6: The transport-mean-free path I* of our model wet 
foam, as a function of the liquid volume fraction e. I* is 
measured in units of the bubble radius R. rn n = 1.0 and 
n ut = 1-34 are assumed. Our analytical result agrees well 
with the relation I* / R ~ 0.11/e + 2.37. Numerical simulations 
of Ref. HH are well described by l*/R « 0.26/e + 4.90. 



some light on the empirical law of Vera et al. and the 
numerical simulation of Sun and Hutzler: l*/R is a lin- 
ear function of 1/e, see Eqs. (HJ and (|23|) . For a better 
understanding of the empirical law |[T[). we aim at a more 
realistic model which not only considers Fresnel's inten- 
sity reflectance with its significant dependence on the in- 
cidence angle, but also the broad distribution of photons' 
step length. Also an extension to the three-dimensional 
packing of polydisperse spheres is envisaged. 
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APPENDIX A: THE PROBABILITY 
DISTRIBUTIONS F(j) AND F'(j') 



immersed in the water (n ou t = 1-34), we find r ~ 0.09 
and I* ~ 9R. Our transport-mean-free paths are smaller 
than the experimental values @, Q by a factor of about 
2. Now we consider a simple model for wet foams. For 
the air bubbles (n,- n = 1) immersed in the water (n ou t = 
1.34), we estimate r w 0.20. Figure [S] delineates I* (in 
units of R) as a function of the liquid volume fraction e = 
1 — 0. From Fig. [5] we find that in the range 0.08 < e < 
0.15, our analytical result agrees well with the relation 



I* 



R(— + 2.37). 



(23) 



Using the hybrid lattice gas model for two-dimensional 
foams and Fresnel's intensity reflectance, Sun and Hut- 
zler performed numerical simulation of photon transport 
[24j . For 0.02 < e < 0.16, their numerical results can 
be fitted to I* w R(0.26/e + 4.90). Fig. [6] compares our 
analytic prediction with the numerical result of Ref. [24| . 
Again our analytic estimate is smaller than the numerical 
simulations by a factor of about 2. 

Quite remarkably, our analytic estimate of the 
transport-mean-free path I* quoted in Eq. ([2|), sheds 



To find the probability distribution of the random 
variable 7, we assume that the impact parameter s in 
Fig. He) bas a uniform distribution in the interval [0, R]. 
In other words, we assume that the number of incident 
rays with an impact parameter less than s is propor- 
tional to s. The cumulative distribution function F c (pj) = 
O F(ip)dip is then F c (-f) = Prob(s < i?sin7) = sin 7. It 
follows that 



F(j) 



c?7 



cos 7. 



(Al) 



Now we consider path of photons inside the disk, see 
Fig. QJd) ■ Each ray can be characterized by its distance 
s' from the center of the disk. We assume that the ran- 
dom variable s' has a uniform distribution in the interval 
[0, R\. Since s' = Rsmj', the cumulative distribution 
function is F' c (j') — sin 7', and 



b (7 ) = — — — = cos 7 . 



(A2) 



Further numerical simulations of our model confirm these 
analytical results for ^(7) and F'(^'). 
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